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Abstract 


We  report  here  some  new  results  we  have  obtained  for  the 
diagonal  element  of  the  second-order  density  matrix  (pair  function) 
and  for  the  momentum  distribution  of  an  interacting  electron  gas. 

In  the  low  density  limit,  the  pair  function  has  been  discussed 
previously  by  March  and  Young  (1959)  and  it  is  shown  here  that, 
to  a  similar  degree  of  approximation,  the  first-order  density 
matrix  and  hence  the  momentum  distribution  may  be  obtained.  In 
the  high-density  problem,  the  pair  function,  has  also  been 
calculated  using  the  Green's  function  approach  and  this  result 
complements  the  calculation  of  the  momentum  distribution  carried 
out  recently  by  Daniel  and  Vosko  (i960).  Prom  the  information 
thus  available  in  these  high  and  low  density  limits,  the  range 
of  usefulness  of  the  concept  of  the  Permi  surface  in  an 
interacting  electron  gas  is  discussed. 

Alsc,  results  of  variational  calculations  on  the  stability 
of  Overhauser  spin-density  waves  in  metals  are  reported .  It  is 
concluded  that  such  waves  reduce  the  energy  below  the  Hartree- 
Pock  value  only  when  the  interaction  becomes  very  strong,  or 
equivalently  when  the  density  is  very  low. 

Some  preliminary  work  on  the  correlation  problem  in  a 
non-uniform  electron  gas  is  also  briefly  discussed. 
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1 ,  Introduction 

As  we  discussed  in  ASR  No.  1,  the  problem  Of  dealing  with 
the  Coulomb  interactions  between  conduction  'electrons  in  metals 
has  provod  very  formidable  and  no  complete  solution  has,  as 
yet,  been  found.  In  recent  years,  beginning  with  the  important 
investigation  of  Macke  (1950),  attention  has  been  focused 
primarily  on  the  use  of  perturbation  theory.  Unfortunately,  the 
perturbation  expansion  does  not  prove  useful  in  practice  in  the 
range  of  real  metallic  densities,  which  are  characterized  by  the 
usual  gas  parameter  r„  lying  in  the  range  2  <  r„  <  5.5. 

As  we  emphasized  in  ASR  No.  1,  in  attempting  to  understand 
the  correlation  problem,  it  is  often  useful  to  consider  a  wider 
range  of  densities  than  that  represented  above,  and  in  particular 
to  consider  the  formal  limiting  cases  r_  -»  0  and  r_  -♦  <».  The 
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first  case  represents  the  domain  of  perturbation  theory  baaed 
on  plane  waves,  but  the  second  case  has  been  less  fully  discussed. 
For  this  reason,  we  have  extended  the  earlier  work  by  March  and 
Young  (1959)  on  the  diagonal  element  of  the  second-order  matrix 
(pair  function)  and  have  thereby  obtained  the  first-order  density 
matrix  and  momentum  distribution  in  this  limit  of  low  density. 

The  results  of  this  calculation  are  reported  in  g3  after  a  brief 
summary  of  the  results  of  the  Hartree-Pock  approximation  in  §2. 

In  §4,  the  other  extreme  limit  of  high  density  is  discussed.  In 
order  to  facilitate  the  calculations  in  this  limit,  a  generalization 
of  the  concept  of  density  matrices  proves  desirable  as  we 
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emphasized  in  §2  of  ASR  No.  1.  The  appropriate  tools  are  then 
found  to  he  the  first  and.  second-order  G-reen’s  functions 
familiar  to  field  theorists.  A  calculation  is  reported  which 
yields  the  pair  function  in  the  high-density  limit;  uniform 
application  of  .  this  same  method  leads,  in  addition,  to  results 
in  agreement  with  those  of  Daniel  and  Vosko  (i960)  for  the 
momentum  distrihution.  By  analysis  of  these  limiting  forms 
for  high  and  low  density  it  is  shown  in  §5  that  information  may 
he  obtained  on  the  range  of  usefulness  of  the  concept  of  a  Permi 
surface  in  metals. 

In  §6  we  report  the  final  results  of  a  variational 
calculation  carried  out  hy  one  of  us  (N.H.M. )  in  conjunction 
with  M.  I.  Darhy  and  W.  H.  Young.  This  was  designed  to  examine 
the  stability  of  Overhauser  spin-density  waves  in  metals.  The 
results  provides  strong  evidence  that  such  waves  are  only 
energetically  favoured  for  very  strong  interactions,  or 
equivalently,  for  very  low  densities.  They  appear  therefore 
to  have  no  relevance  to  the  correlation  problem  for  real  metallic 
densities. 

Finally,  in  §7,  we  report  some  preliminary  findings  on  the 
proton  lattice,  in  both  high  and  low  density  limits. 

2.  Hartree-Pock  pair  function  and  first-order  density  matrix 

We  briefly  recall  here  that  our  problem  is  basically  to 
calculate  the  gzround  state  wave  function  £2  •  •  •  %)  for 
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a  system  of  N  electrons,  where  denotes  the  space  co-ordinates 
r^  and  the  spin  co-ordinates  from  the  SchrBdinger  equation 

H\|r  =  (2.1) 

where  the  Hamiltonian  H  is  given  in  full  in  equation  (2.1)  of 
ASR  No.1.  As  stressed  in  that  report,  it  is  often  "both  convenient, 
or  more  important,  conceptually  very  fruitful,  to  work  with  first 
and  second-order  splnless  density  matrices  Y(ri )  and 
r(21i '£.2 'rir2 )  defined  by  the  equations 

y(£i '£1  )  =N  |\lf*(ri '0iX2...Xjj)t(£i0i  dX2  . .  .d2Cjj  (2.2) 

r(£l  '£2'£i£2)  =  |t*(£l  '0ir2'02X3...3^)t(£l0l£2O2X3.  . .%) 

dOi  d02  dxa  ...  d^  (2.3) 

We  note  now  that  the  diagonal  elements  of  y  and  r  have  direct 
physical  interpretations.  Thus  y(ri  ri  )  gives  the  particle 
density  which  in  the  Sommerfeld  model  is  simply  the  constant 
3/47cr_®,  while  r(rir2ri£2)  is  effectively  the  probability  of 
electron  separation  or  the  pair  distribution  function.  Then, 
as  Mayer  (1955)  was  the  first  to  show,  the  energy  per  particle 
e/N  is  given  by 

-1  * 

ii’(e'c)]  -  2^  J  (1  -  F(S))5  «  (2.4) 

/ 


where  P(J^)  =  P(lil'  “  2,|)  is  the  pair  fxmction,  normalized  to  unity 
for  large  | r'  -  r| . 

When  the  total  wave  function  is  a  determinant  of  plane  waves, 
we  have  the  usual  Hartree-Pock  forms 


Y(r'r)  = 


k  ® 

Kf 


(k-lr'  -  r  [) 

'’k'^  r~  -rT  p  -  p  cos  p)/p®  (2.5) 


P(r'r)  =  1 


(2.6) 


where  k^  is  the  magnitude  of  the  wave  vector  at  the  Perml  surface 
and  is  related  to  the  mean  interparticle  spacing  by 


‘'f  ■'a 


■w 


(2.7) 


Using  atomic  \inits,  it  follows  from  (2.5)  and  (2.6)  that  the  energy 
per  particle  is 


e 

N 


(2.8) 


3.  Pair  function  and  first-order  matrix  in  low-density  limit 

The  low-density  form  of  the  pair  function  has  been  discussed 
earlier  by  March  and  Young  (1939)  using  the  electron  lattice  model 
of  Wlgner  (1938).  We  simply  summarize  the  results  by  saying  that 
as  r_  becomes  very  large,  the  electrons  relative  to  a  given  particle 

O 

we  have  singled  out,  and  placed  at  the  origin  of  co-ordinates  sit 
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on  the  sites  of  a  body-centred  cubic  lattice.  As  is  then 
reduced  somewhat,  the  electrons  can  be  thought  of  as  represented 
by  harmonic  oscillator  functions 


The  results  obtained  by  March  and  Young  (1959)  for  the  pair  function 
are  represented  in  curves  1  and  2  of  Pig,  1,  for  cases  r_  =  100 
and  r_  =  4.  For  comparison,  the  Fermi  hole  result  given  in  (2.6) 
is  also  shown  in  curve  3*  We  emphasize  that  while  curve  1  should 
be  reliable,  curve  2  represents  an  extrapolation  beyond  the  range 
of  validity  of  the  low  density  form  (3.1),  and  is  given  solely  to 
show  the  qualitative  influence  of  varying  r„. 

We  shall  now  indicate  how  this  work  may  be  generalized  to 
yield  the  first-order  density  matrix  and  momentum  distribution  for 
a  low-density  gas.  Prom  the  orbitals  (3*1),  centred  on  each 
lattice  site,  we  build  a  Dirac  density  matrix.  This  will  not,  of 
course,  contain  the  condition  of  translational  invariance,  that  is 
the  diagonal  element  will  not  be  a  constant.  To  obtain  a  result 
consistent  with  the  Somraerfeld  model  we  must  average  over  all 
positions,  and  then  the  desired  first-order  density  matrix  is 
easily  shown  to  be 

Y(r'r)  =  ^  exp  g|r'  -  r|®^.  (3.2) 

The  form  (3.2),  as  is  easily  seen,  satisfies  all  the  essential 
conditions,  but  its  range  of  validity  is,  of  course,  restricted. 


-6- 


as  we  shall  discuss  helow.  Nevertheless,  In  conjunction  with  the 
high  density  results  of  §U,  it  may  he  used  to  draw  some  interesting 
conclusions  about  the  Perml  surface  (see  §5). 

At  this  stage,  it  is  enlightening  to  examine  the  momentum 
distribution  corresponding  to  the  first-order  matrix  (3.2),  and 
this  may  be  found  as  follows;  We  require  the  occupation  numbers 
P(i)  of  plane  wave  states  e^*^,  where  V  is  the  volinne  of  the 
metal  and  thus  we  write 

y(e'e)  =  I  y  P(k)e-‘^-£'  (3.3) 

k 


or  remembering  that  the  density  of  states  in  k  space  is  (l/8x®)V 
and  using  Bauer's  expansion  for  a  plane  wave  as  a  series  of 
spherical  waves: 

.  sin  kir'  -  r  j 

y(£'e)  «  Jjp-  J  P(k)  dk,  (3.4) 

It  is  convenient  at  this  point  to  measure  k  in  units  of  the  Fermi 
momentum  k^,  that  is  we  write  K  s  and  then  (3*4)  becomes 


sin  k^K)r'  -  £) 


K^dK, 


(3.5) 


Inverting  this  relation  and  using  (3.2)  for  yCe'e)  we  find 

P(K)  =  ^  exp  j^- 


(3.6) 


where  we  have  eliminated  using  (2.7)*  As  we  have  remarked 
earlier,  the  range  of  validity  of  (3*6)  is  restricted,  hecause 
the  orbitals  t  on  different  lattice  sites  hav^  been  assumed 
orthogonal,  whereas  this  Is  only  rigorously  true  in  the  limit 
r  -♦  00.  However,  a  rough  estimate  of  the  range  of  validity  may 
be  obtained  by  noting  that  the  occupation  numbers  P(K)  must  always 
lie  between  0  and  1,  Since  P(K)  as  given  by  (3.6)  has  its 
maximum  at  K  =  0,  we  must  have  that 


(3.7) 


or  rg  >  9*3. 

The  kinetic  energy  per  particle,  in  the  approximation  implied  by 
(3.6)  is  easily  shown  to  be 


T 

N 


(3.8) 


in  agreement  with  Wigner  (1938).  La  er  treatments  based  on  more 
careful  study  of  the  dynamics  of  the  vibrating  electron  lattice 
(Coldwell-Horsfall  and  Maradudin,  1960;  Carr,  1961)  change  (3*8) 
somewhat,  but  will  not  affect  the  overall  validity  of  the  argument 
presented  above.  We  note  finally  for  the  low-density  limit, 
that  the  momentum  distribution  is  of  classical  'Maxwell-Boltzmann* 
form:  this  is  in  sharp  contrast  to  the  high  density  results  to 


which  we  now  t\irn 
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4.  Pair  function  and  first-order  density  matrix  In  hlgh-denslty 
limit 

Using  the  perturbation’ theory  of  Gell-Mann  and  Bnjeckner 
(1957)  in  suitable  form,  we  shall  now  discuss  the  way  in  which  the 
pair  ftinction  and  the  first-order  density  matrix  develop  from  the 
Hartree-Pock  forms  given  in  ^2,  equations  (2.6)  and  (2.5).  In 
order  to  do  so  it  is  very  convenient  to  introduce  a  generalization 
of  density  matrices  to  include  also  the  time;  •  this  will  greatly 
facilitate  the  setting  up  of  useful  equations  from  which  the 
matrices  can  be  determined.  We  discussed  this  situation  fully 
in  §2 (ill)  of  ASR  No.  1,  where  we  pointed  out  that  the  first  and 
second  order  Green's  fiinctlons  constitute  the  desired  generalizations 
of  the  density  matrices.  In  fact,  equations  (2.38)  and  (2.39)  of 
ASR  No.  1  give  the  explicit  relations  between  the  Green's  functions 
and  the  density  matrices.  The  basic  equation  for  the  second-order 
Green's  function  (see  (2.35)  of  ASR  No.  1)  may  be  written  in  the  form 

0(XiXaX3X4)  =  G(XiX3)G(X2X4)  -  G(xiX4)G(X2X3) 

+  j  d^yi  d^^ya  d'^Zi  d^^Za  G(Xiy,  )G(xaya)l(yiy2ZiZa)0(ziZaX3Xi). 

(4.1) 

In  the  high  density  limit  it  now  appears  possible  to  make  progress 
in  the  electron  gas  problem  by  Inserting  the  first  approximation 
for  G(xiXaX3X4)  Inside  the  integral,  provided  we  approximate 
sufficiently  carefully  to  the  interaction  operator.  The  procedure 
to  be  adopted  was  outlined  in  §4(iil)  of  ASR  No.  1,  and  we  have 
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now  carried  through  the  calculations  in  detail,  as  described  below. 

The  basic  equations  of  the  treatment  we  shall  employ  are  (4.^5)  and 
(4.26)  of  ASR  No.  1,  The  method  used  is  to  seek  the  Fourier 
transform  of  the  two-particle  Green's  functioil,  which  is  closely  related 
to  the,  second  order  density  matrix  (see  (2,39)  of  ASR  No.  1), 

We  therefore  write  the  interaction 

v(r)  =  ^  =  i  v(k)  e^^*^  (4.2) 

and  then  it  follows  that 

v(k)  =  (4.3) 

Furthermore  we  develop  Go(rt),  which  enters  the  Green's  function 
equation  for  the  two-spin  case,  namely, 

G(rt£' trt"r' t")  =  Go(rt  rt")Go(r't£'t")  -  ^  Go(£t£'t" )&©(£' trt") 

-  J  d'^yi  d^ya  Go(Etyi  )Go(ji'tya)V(yi-ya)Go(yirt")Go(yaE't") 

+  i  J  d^ya  Go(rtyi  )Go(r'  ty2)v(yi-ya)Go(yi£'t")Go(y2£t"), 

(4.4) 

in  the  form 

00 

Go(rt)  =  ^  y  I  Go(ke)e^(^**’  "  dS 

k 

=  7  y'  Oo(kt)e^^‘^ 


(4.5) 
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After  some  calculation  it  may  then  ho  ahovm  that 

Go(kt)  =  -  * 

=  0 


t  <  0  k  >  0 

t  >  0  k  >  1 

t  <  0  k  >  1 


(4.6) 


t  >  0  k  <  1 

Also  it  is  known  that  equation  (4.25)  of  ASR  No.  1  for  the 
modified  interaction  V  may  be  solved  by  Fourier  analysis,  and 
writing 

V(yi-y2)  =  2^  J  v(q<D)  da>  (4.7) 


the  solution  is 


v(qa))  = 


2TOcr< 


Qo(qw)3 


(4.8) 


where 


Q,(»«)  =  [  d>k  L  ■  I  a.lt)  ■  »  +  (aVsV  q.k)+I*j 

■kl<1  '• 

k+ql>1  (4,9) 

This  latter  function  Qo(q,w)  has  been  discussed  earlier  by  DuBois 
(1959). 

If  we  now  denote  the  correction  to  the  Hartree-Pock  result  for 
0(rt  r't  rt"  rt")  by  Gj^g,  then  we  may  write 
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Gjjg(rt  r't,  rt"  r't")  3  jdq. 

=  J  d'^yi  d^^ya  Go(rt,  y^ )  Oo(r't  ya)  V(yi-ya)  Go(rirt")Go(rar' t") 
"  'fa'kV’*^  /  ^<1  ^  gigCr-r')  -  i(e+e')(t'-  t") 


X  Go(ke)  Go(k  -  g,  e  -  a>)  v(qco) 

If  we  let  t"  -♦  t"*"  at  this  stage  we  obtain 


(U.10) 


dk  d€  dk'  de'  dco  Go(ke)  Go(k  - 


<lf 


6—0)) 


v(qo))  Go(k'e')  Go(k'  +  g,  6'  +  co) 


=  -  ^2^)7  J  Qo^(gw)  v(go))  do>.  (4.11 ) 

Substituting  for  v(go))  we  obtain  the  final  result 

-  /  N  1  f*  Qo^(go))  arsA°g°  . 
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Further  progress  calls  for  the  use  of  numerical  methods  and  this 
aspect  of  the  work  is  now  being  planned. 


5.  The  Meaning  of  a  Fermi  sxirface 

The  fact  that  we  can  obtain  the  momentum  distribution  in  two 
limiting  cases  has  interesting  implleatlons.  For,  as  the  work 
of  Daniel  and  Vosko  has  shown,  the  discontinuity  In  P(K)  at  the 


-12 


Fermi  surface  for  r  =  0  is  not  removed  but  only  reduced  as  the 

O 

density  is  lowered.  While  this  result  seems  not  to  have  been 
established  rigorously  outside  perturbation  theory,  there  seems 
strong  empirical  evidence  from  the  sharpness  of  the  Fermi  surface 
in  real  metals  that  such  a  discontinuity  does  in  fact  persist. 
However,  our  low  density  treatment  r_  oo  could  equally  well  be 

O 

made  the  basis  of  a  perturbation  treatment,  and  no  sign  of  the 
discontinuity  would  then  occur.  Thus,  the  evidence  seems  clear 
that  as  we  follow  the  momentum  distribution  as  it  develops  from 
the  two  limiting  cases  of  small  and  large  r_  there  must  come  a 

O 

critical  coupling  strength,  or  a  critical  density,  at  which  the 

discontinuity  in  the  momentum  distribution  is  reduced  to  zero. 

For  lower  densities,  it  then  appears  that  the  concept  of  a  Fermi 

surface  will  no  longer  be  useful.  Vi/e  believe  at  present  that  the 

critical  value,  r  say,  will  lie  outside  the  range  of  real  metallic 

c 

densities,  or  in  other  words  that  r.  will  exceed  5*5»  No 

c 

quantitative  evaluation  has  so  far  proved  possible  however. 
Questions  also  remain  as  to  the  nature  and  order  of  the  'transition' 
occurring  at  r_.  We  do  not  expect  that  the  pair  function  will 
undergo  any  marked  changes  at  r.  and  probably,  to  obtain  a  clear 
picture  of  the  transition,  it  will  eventually  prove  necessary  to 
examine  the  third  and  higher-order  density  matrices. 

Finally,  to  give  an  indication  of  the  change  in  the  momentum 
distribution  as  the  density  is  varied,  we  have  plotted  in  curve  1 
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of  Figure  2  results  for  P(k)  as  given  by  Daniel  and  Vosko  for  the 
high  density  values  (rg  =  2)  and  by  (3.6)  for  the  low  density 
region.  Curve  2  for  r^  =  16  should  be  viewed  with  some  caution, 
but  the  result  for  r_  =  100  (curve  3)  should  represent  a  good 
approximation.  The  groat  difficulty,  as  we  have  stressed,  is  to 
obtain  a  sufficiently  good  approximation  to  the  ground  state  in 
the  region  2  <  r_  <  5*5  and  to  locate  the  critical  value  of  r_ 
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with  precision. 

6.  Stability  of  spin  density  waves  in  metals  (with  M.  I.  Darby  and 
W,  H.  Young) 

Wo  record  here  the  final  results  of  the  calculations  of 
M.  I.  Darby,  one  of  us  (N.H.M.)  and  W.  H.  Young  bearing  on  the 
spin  density  waves  of  Overhauser  (1960).  The  analysis  was  set 
out  in  ASR  No.  1  and  more  accurately  (avoiding  a  sphere 
approximation  made  earlier)  in  section  2  of  TSR  No.  4.  However, 
evaluation  of  the  potential  energy  given  in  equation  (3)  of  that 
Report  indicated  that  the  term  in  a*  in  the  potential  energy  was 
exceedingly  small.  Examination  of  the  reasons  for  this  led  to 
the  use  of  an  alternative  variational  function  which  proved  much 
more  accurate. 

Using  the  notation  of  ASR  No.  1,  the  transformation  function 
E(£)  taken  as  having  components 

Rj  =  Xi  +  I  Sin  a  ^')  (6.1) 
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and  then  in  place  of  equation  (1)  of  TSR  No,  4»  we  easily  find 
for  the  increase  in  the  kinetic  energy  per  particle: 


Similarly  the  change  in  the  potential  energy  per  particle  is 
given  hy 


(6.3) 


Thus,  the  energy  is  lowered  by  allowing  non-zero  amplitude  of 
the  spin-density  waves  when 

X  >  187 

or  rg  =  X  >  116.  (6.4) 


This  appears  to  us  to  be  fairly  conclusive  evidence  that  spin 
density  waves  are  not  energetically  favovirable  until  the  electron 
density  becomes  very  much  lower  than  occurs  in  real  metals.  Charge 
density  waves  may  perhaps  lower  the  energy  at  higher  densities 
however,  (see,  for  example,  Yoving,  1961). 


7.  Proton  lattlco  problem  (with  W,  Jones) 

We  shall  now  give  a  brief  accoxint  of  some  preliminary  work 
we  have  carried  out  on  the  non-uniform  elcctijon  gas  problem, 
'JVhereas,  in  the  Sommerfeld  model,  the  'natural'  orbitals,  of  the 
theory  are  plane  waves  for  all  electron  densities,  in  the  case  of 
a  periodic  lattice,  the  'natural'  orbitals  will  vary  with  density 
This  is  an  added  complication,  and  has  led  us  to  investigate  how 
we  should  obtain  the  Bloch  orbitals  for  a  given  density.  This 
seems  possible  at  present  only  in  the  limits  of  high  and  low 
density,  and  we  shall  restrict  our  discussion  to  these  cases. 

We  should  remark  at  this  point  that  the  ground-state  energy  of 
an  electron  gas  in  a  lattice  of  positive  point  charges  has  been 
studied  by  Bellomans  and  de  Leaner  (196I).  These  workers 
expanded  the  grand  partition  function  of  the  system  in  powers  of 
the  coupling  parameter  between  the  electrons  and  the  positive 
point  charges  and  between  pairs  of  electrons.  By  summing 
infinite  series  of  divergent  terms  they  were  able  to  obtain 
finite  results  for  the  energy  per  particle. 

We  have  restricted  ourselves  up  to  the  present  to  the 
calculation  of  Dirac's  matrix  in  the  high-density  limit.  Thus, 
we  cannot  expect  to  obtain  an  energy  which  includes  the  Gell-Mann 
Brueckner  correlation  energy,  in  contrast  to  Bellemans  and 
de  Leener,  However,  we  shall  demonstrate  that  there  is  an 
Intimate  connection  between  their  approach  and  that  given  below. 
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We  shall  work  in  a  Hartreo  approximation  at  first,  and  shall 
write  the  periodic  self-consistent  potential  V(r)  in  the  form 


V(r)  =  y  V  e 

L.  ^ 

^£n 


“ 


(7.1) 


where  the  ^  are  reciprocal  lattice  vectors.  We  now  make 
use  of  the  first-order  equation  of  March  and  Murray  (1961),  which 
we  may  write  as 


k-®  k_®  ^  r 

E  \J 


ISn-E  „  f  ISn-E 


r —  dr  (7.2) 


where  the  Fermi  wave  vector  has  magnitude  k^  and  is  related  to 
the  Fermi  energy  ^  hy  k^*  =  2^.  Now  if  the  parameter  \  measures 
the  magnitude  of  the  positive  point  charges  at  the  lattice  sites 


R„ ,  then 
-n’ 


7=V(r)  =  X  d(r  -  R^)  - 


4np 


Sn 


-  Uicp, 


(7.3) 


^n 


where  H  is  the  volume  of  a  unit  cell.  Hence,  substituting  (7.1) 
and  (7.2)  in  (7.3)  we  may  write 
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-E 


i/c_  .r 


/c  ®  V  e 


-  4« 


( 


k  3 
3^ 


Ir  2 

2t^ 


l-  \h 


i^Sy^.H  di(2k^r) 


r —  dr 


■r). 


(7.4) 


Hence,  for  ^  ^  0,  we  have 
—n  ' 


.  r  _  ,  2k.=  !  i£n-E  J'tSkfl-)  _  ^  - 

!%  i  =  ■  n- 


Now 


/ 


:Ji(2k-r) 

3^  =(^j  /  a 


•  ■  W*  /  - 


T 


^  4ii^  / 


2k^r  -  2k^r  cos  2k^r} 


dr 


1  n 


(7.6) 


where 


J(kf ,  Kj,)  =  +  (kf.=*  -  -g-) 


■ftn 


2k^  - 


(7.7) 


To  begin  our  discussion  of  the  total  energy,  we  first 
investigate  the  E(k)  relation  in  the  present  approximation. 
Thus,  we  write 


E(is)  =■  y 


”11 


(7.8) 
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where 


V  =  /  Hp(r'  -  r)  dri, 


(7.9) 


with  the  rule  that  the  Hamiltonian  H  operates  only  on  r  and  we 
put  r'  =  r  before  performing  the  integration.  p(r',  £)  here 
is  the  density  matrix  for  the  filled  band.  After  some 
calculation  we  then  find,  writing 


E(k)  =  ^  +  AE(k), 


(7.10) 


that 


(7.11) 


‘m 


where  is  the  volume  of  the  Brillouin  zone  and  kjjj  = 

To  obtain  the  change  in  the  sum  of  the  eigenvalues  we  rogpiire 

Ae  =  ^  f  AE(k)d^ 

^  ^  Fermi  sphere 


m 


■m 


Fermi  sphere 

(7.12) 

The  integral  may  be  evaluated  in  terras  of  the  fxinction  J(k^,  /Cj^^) 
defined  earlier,  and  the  result  may  be  written 


(7.13) 


m 
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If  wc  now  insert  the  Fourier  components  of  the  unscreened  lattice 


V. 


we  find,  introducing  u  = 


Ae 


is  wo  put 

it 

CM 

-  u®)  , 

1  +  u[ 

1  -  u| 

(7.14) 


(7.15) 


u 


This,  in  fact,  agrees  with  the  last  term  in  the  energy  as  given 
hy  Bellemans  and  de  Leener,  and  means  that  their  energy  can  he 
written  down  directly  by  straightforward  methods,  without  recourse 
to  their  development  of  the  grand  partition  function.  Obviously, 
our  result  (7.13)  goes  beyond  their  approximation  if  we  insert 
the  correct  Fourier  components. 

We  finally  note  that  in  the  low  density  limit,  the  hydrogen 
atom  wave  functions  must  give  ^he  solution  of  the  problem.  Then 
it  follows  that  the  momentum  distribution  function,  which  seems 
of  some  importance  in  following  the  behaviour  of  the  system  as 
the  lattice  parameter  is  varied,  has  the  form 


It  would  certainly  be  of  interest  to  enquire  what  modifications 
would  occ\ir  in  this  function  as  the  lattice  parameter  is  now 
decreaded,.  Obviously,  a  better  approximation  than  e”**  to  the 
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localized  orbitals  is  then  required,  and  we  are  at  present 
considering  ways  of  introducing  perturbations  having  the  correct 
crystal  symmetry  in  order  to  obtain  the  form  of  l(p)  away  from 
the- point  at  infinity. 

8.  Conclusion 

For  the  Somraerfcld  model,  the  physical  behaviour  of  the 

system  in  the  two  limiting  cases  r^  -»  0  and  -►  «  Is  now  well 

established.  If,  as  we  believe,  a  transition  occurs  at  a 

critical  density  r  ,  then  accurate  location  and  description  of 

c 

this  is  obviously  of  considerable  Interest  and  we  are  continuing 
to  examine  this  problem. 

Less  is  known  on  the  non-tiniform  gas  problem  particularly 
in  the  low-density  limit,  and,  as  remarked  above,  we  are 
investigating  crystalline  field  methods  of  developing  (7.16) 
for  higher  densities. 
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Captions  for  Flguroe 


Page 

Pair  functions  for  an  electron  gqs.  2i 

Curve  1.  Low  density  form  for  r^  =  100. 

Curve  2.  Extrapolated  low  density  form  for  r_  s  U. 

O 

Curve  3.  Hartree-Pock  pair  function  of  eqn.  (2,6), 
correct  in  limit  r^  -♦  0, 


Momentum  distribution  in  an  electron  gas. 

Curve  1 .  Daniel  and  Vosko  results  for  r^  =  2, 

Curve  2.  Low  density  form  for  r^  =  16. 

Curve  3,  Low  density  form  for  r_  =  100, 
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